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ABSTRACT 

High  Frequency  Surface  Wave  Radar  ( HFSWR )  takes  advantage  of  the  diffraction  of  electromagnetic 
waves  over  the  conducting  ocean  surface  to  detect ,  locate,  and  track  surface  vessels  beyond  the  line- of  - 
sight  horizon  of  the  earth.  However,  long-range  surface  vessel  detection  is  often  confounded  by  radar 
clutter  comprising  echoes  from  the  ionospheric  plasma.  In  this  paper,  we  characterize  the  angular 
spectrum  of  these  reflections,  and  from  this  information  deduce  the  signal-to- clutter  processing  gain 
that  can  be  realized  by  various  adaptive  receive  antenna  array  configurations.  In  particular,  it  is  shown 
that  there  is  advantage  to  be  realized  by  sampling  ionospheric  echoes  with  a  planar  two-dimensional 
array  rather  than  a  conventional  linear  one- dimensional  HFSWR  array.  Using  a  planar  array,  the 
radar  can  distinguish  high- elevation  ionospheric  clutter  signals  from  low-elevation  surface  target  echoes. 


1  ADAPTIVE  ARRAYS 


Ionospheric  reflection  in  HFSWR  applications  complicates  the  detection  of  weak  target  echoes.  In 
particular,  during  nighttime  Spread  F  conditions,  ionospheric  clutter  can  occupy  most  of  the  far 
ranges  in  the  radar  range-Doppler  search  plane.  Adaptive  antenna  array  techniques  can  be  used  to 
increase  the  ratio  of  surface  vessel  echo  power  to  ionospheric  clutter  power  in  cluttered  radar  ranges. 

We  consider  a  surface  vessel  radar  echo  consisting  of  a  plane  wave  of  amplitude  x  that  is  received 
by  an  array  of  antenna  elements  arranged  in  an  arbitrary  configuration.  The  amplitude  and  phase 
response  of  the  individual  sensors  to  the  plane  wave  is  captured  by  an  array  manifold  vector  v.  In 
addition  to  the  plane  wave  signal,  the  sensors  are  corrupted  by  a  “noise”  vector  n  that  is  dominated 
by  ionospheric  clutter  signals.  Thus  the  signal  measured  by  the  array  can  be  written  as  y  =  vx  +  n. 


We  assume  the  radar  signals  are  narrowband  such  that  they  can  be  described  by  complex  amplitudes. 
We  devise  a  linear  estimator  for  x  that  consists  of  summing  the  measured  signals  according  to  a 
“weight”  function  w.  This  estimator  is  denoted  as  x  =  w^y,  where  H  is  the  complex  conjugate 
transpose.  In  order  to  minimize  the  variance  of  the  expected  value  E{x)  under  the  constraint  of  no 
distortion  E(x)  —  E(x),  the  optimal  choice  of  weights  in  this  linear  estimator  is  given  by  [1]: 


w 


H 


v^R-1 

v  ±v"n 

v^RnV 


(1) 


where  Rn  =  £(nn^)  is  the  noise  covariance  matrix.  This  estimator  is  an  adaptive  beamformer,  pro¬ 
viding  constructive  interference  in  the  target  direction  and  destructive  interference  in  noise  directions. 
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ORGANIZATION 


This  beamformer  is  distortionless,  so  it  provides  no  improvement  in  target  signal  power.  All  SNR 
improvement  arises  from  reducing  noise  power  associated  with  the  ionospheric  clutter.  At  the  input 
to  the  beamformer,  the  noise  power  received  by  individual  antenna  elements  is  given  by  the  diagonal 
entries  of  Rn.  To  help  quantify  the  reduction  in  this  noise  power  by  the  beamformer,  we  assume  the 
noise  power  is  the  same  at  all  sensors,  and  equal  to  a  scalar  denoted  by  Rn.  The  noise  power  at  the 
output  of  the  beamformer  is  given  by 


Rj.  —  w^Rnw 


1 

v*R^V 


(2) 


Therefore,  the  SNR  gain  due  to  the  array,  denoted  as  the  “array  gain” ,  is  given  by 


A  =  Rn/Rx  =  Rn^YiJw  =  vHpn1v. 


(3) 


Here,  pn  is  a  normalized  covariance  matrix  where  the  entries  on  the  matrix  diagonal  are  unity.  In  the 
case  of  uncorrelated  noise,  the  array  gain  will  be  simply  A  =  vHv  =  N ,  where  N  is  the  number  of 
sensors.  In  the  case  of  correlated  noise,  the  array  gain  could  be  more  or  less  than  N. 


2  DIRECTIONAL  SPECTRUM  OF  IONOSPHERIC  ECHOES 


To  estimate  pn  for  arbitrary  receive  arrays,  we  can  characterize  the  ionospheric  echo  as  a  spectrum 
of  plane  waves  5(k,a;)  incident  on  the  array,  where  k  is  a  vector  wavenumber  (kx,ky,kz),  and  uj 
is  the  wave  frequency.  In  the  case  of  free  space,  *S(k,  uj)  reduces  to  a  two-dimensional  distribution 
of  horizontal  wavenumbers  S{kx,ky )  since  uj  is  fixed  at  the  radar  carrier  frequency  and  the  vertical 
wavenumber  kz  is  equal  to  kz  —  ^ ( uj/c )2  —  k%  —  k%  by  the  freespace  dispersion  relation  uj  =  c|k|. 
S{kx ,  ky )  is  equivalent  to  the  distribution  of  wave  direction  of  arrival  (DOA)  in  elevation  and  azimuth. 


To  estimate  S{kx ,  ky )  we  extend  to  two  dimensions  a  method  that  was  devised  for  the  characterization 
of  one-dimensional  plasma  turbulence  using  fixed  probe  pairs  [2].  The  plane  wave  spectrum  S(kx ,  ky ) 
is  related  to  the  plane  wave  spatial  autocorrelation  function  R  by 

Jl  dkx  dky  S(kx,  ky)eik*x+ikyY  =  R(X,  Y),  (4) 

where  X  and  Y  are  the  lags  in  the  x  and  y  directions,  and  the  autocorrelation  function  R  is  defined 
as 

R(X,Y)  =  E[**(x,y)$(x  +  X,y  +  Y)].  (5) 

Here,  $  is  a  wave  field  with  random  amplitude  a  and  random  phase  9: 

$(x,y)  =  a(x,y)eld(x’v\  (6) 


In  analogy  to  instantaneous  frequency  in  the  time  domain,  we  define  local  wavenumbers  Kx(x,y)  and 
Ky(x,y)  in  the  spatial  domain: 


Kx(x,y) 


dOjxj,  y) 

dx 


Ky{x,y) 


dd{x, y) 
dy 


By  the  moment-generating  property  of  the  Fourier  transform,  we  have  that 


1  rr  dMdN 

^2  JJ  dkx  dky(ikx)M (iky)N S(kx,ky)  =  ■^W^WR{X,Y) 


(27 r 


X,Y= 0 


(7) 

(8) 
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=  E 


a(X,y)e-^ZE^MVN]  ^ 


m= 0 n=0 


m\n  dXmdYn 


cl{x+X,  y+Y)  ev 


-I  X,Y= 0 


If  the  amplitude  and  wavenumber  spatial  variation  is  slow  compared  to  a  wavelength  such  that 


dma(xi  y ) 


a(x,y)  dXm 


1 


d  mKx(x,y) 


Kx(x,y)  dX ™ 


«  \K?(x,y)\ 


1  dna(x,y) 

1  dnKy(x,y) 

Kn(-r 

a(x,y)  dYn 

•> 

Ky{x,  y)  dY - 

y  V ) 

for  all  ra,  n  >  0,  then  it  can  be  shown  that  Equations  (8)  and  (9)  reduce  to 

^2  JJ  dkxdkyk^k^S(kx,ky)  «  E[a\x,y)K^(x,y)K^(x,y)]. 


(27 r) 


(9) 

(10) 

(11) 

(12) 


Since  a,  Kx ,  and  are  random  variables,  they  have  a  joint  probability  density  function  (PDF) 
f(a2,Kx,Ky ): 


(27 r) 


JJ  dkxdkyk^ky S(kx,ky)  «  JJJ  da2  dKx  dKy  a2(x,y)K^ (x,y)Ky  (x,y)f(a2 ,  Kx,  Ky),  (13) 


and  so  the  plane  wave  spectrum  in  the  vicinity  of  the  receive  array  can  be  written 

S(fcx,  ky)  «  (2tt)2  J  da2  a2(x ,  y)f{a2,  Kx,  Ky).  (14) 

Thus  the  wave  spectrum  is  proportional  to  a  power- weighted  joint  PDF  of  the  local  wavenumbers.  This 
result  is  analogous  to  the  time  domain  result  that  the  frequency  spectrum  is  equal  to  the  spectrum  of 
instantaneous  frequencies  when  the  signal  modulation  rate  is  slow  compared  to  the  signal  bandwidth. 


3  DIRECTIONAL  SPECTRUM  MEASUREMENTS 

In  practice,  the  right  side  of  Equation  (14)  can  be  estimated  by  compiling  a  power- weighted  histogram 
of  the  phase  lags  from  antenna  pairs  spaced  in  the  x  and  y  directions.  To  this  end,  experiments 
were  carried  out  on  15  Oct  2003  at  an  HFSWR  installation  at  Cape  Race,  Newfoundland,  Canada. 
This  facility  is  described  in  general  terms  in  [3].  For  the  characterization  of  the  ionospheric  signals, 
the  facility  was  operated  in  a  ship  surveillance  mode  at  3.2  MHz,  and  ionospheric  Spread  F  clutter 
signals  were  collected  by  monopole  receive  antenna  pairs  that  were  spaced  54.0  m  in  the  x  direction 
(approximately  southeast)  and  33.3  m  in  the  y  direction  (approximately  northeast).  The  data  were 
analyzed  for  phase  lags  and  power  content  and  the  results  are  plotted  in  Figure  1. 

The  left  panel  of  Figure  1  shows  a  set  of  DO  A  measurements  where  zenith  is  the  origin  of  the  plot  and 
the  horizon  is  the  edge  of  the  circle.  Here,  Kx  is  the  horizontal  axis,  and  Ky  is  the  vertical  axis.  This 
scatter  plot  suggests  that  S(kx,ky )  depends  on  the  radial  distance  yTc |  +  rather  than  kx  and  ky 
individually.  The  center  and  right  panels  are  plots  of  the  marginal  distributions  S(kx )  =  /  dky  S(kx ,  ky) 
and  S{ky)  =  f  dkx  S(kx,ky),  which  take  the  form  of  Gaussian  distributions  of  comparable  variance. 
These  experimental  results  suggest  that  we  model  the  power  spectrum  S{kx,ky )  as  an  uncorrelated 
bivariate  Gaussian  distribution  centered  at  approximately  zenith: 

27T  ^+ky 

S(kx,ky)  «  .  (15) 
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Figure  1:  Left,  measured  DOA  scatter  plot.  Center,  measured  S(kx ).  Right,  measured  S(ky). 

We  will  retain  the  standard  deviation  a  of  this  distribution  as  a  parameter  in  our  analysis.  The 
parameter  cr  is  a  function  of  the  roughness  of  the  reflecting  ionospheric  plasma,  and  describes  the 
angular  spread  a  of  the  ionospheric  reflections: 


<a  =  2sin  1(a/k).  (16) 

The  plots  of  Figure  1  suggest  an  angular  spread  a  of  about  16  degrees. 

To  compute  the  normalized  array  covariance  matrix  pn,  we  evaluate  the  spatial  autocorrelation  func¬ 
tion  R  associated  with  the  plane  wave  spectrum  S(kx,  ky)  at  points  corresponding  to  element  locations 
in  the  receive  array.  The  spatial  autocorrelation  function  is  given  by 

R(X,  Y )  =  ^2  /  dk -  dky  S KYkxX+ikvY  =  (17) 

Thus  any  pair  of  array  elements  in  a  two-dimensional  array  have  a  correlation  that  decays  as  a  Gaussian 
with  respect  to  their  separation  in  the  x  and  y  directions,  where  the  correlation  distance  is  1/a.  In 
this  experiment  the  correlation  distance  is  about  100  m. 

4  ADAPTIVE  ARRAY  GAIN  CALCULATIONS 

We  can  now  determine  the  array  covariance  matrix  pn  for  various  array  configurations  by  evaluating 
the  spatial  autocorrelation  function  R  at  the  locations  of  the  array  elements.  The  adaptive  array  gain 
A  for  these  array  configurations  can  then  be  calculated  according  to  Equation  (3).  This  procedure 
has  been  carried  out  for  four  common  array  configurations,  each  of  N  =  16  elements  and  interelement 
spacing  of  A/2.  These  configurations  consist  of  (1)  a  16-element  linear  array,  (2)  a  cross  of  two  8- 
element  linear  arrays,  (3)  a  16-element  circular  array,  and  (4)  a  filled-aperture  4x4  square  array.  The 
variation  in  A  has  been  calculated  as  a  function  of  the  azimuthal  angle  of  the  target  signal  with  respect 
to  array  broadside,  and  as  a  function  of  the  angular  spread  a  of  the  ionosphere  signals,  as  given  by 
Equation  (16).  The  results  are  plotted  in  Figure  2. 
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Figure  2:  Left,  adaptive  array  gain  versus  target  azimuth  angle  with  respect  to  broadside  (clutter  spread 
angle  =  16  degrees).  Right,  adaptive  array  gain  for  broadside  target  versus  clutter  spread  angle. 

The  results  have  physical  interpretations  as  follows.  The  adaptive  array  gain  for  a  fixed  clutter  spread 
angle  (left  panel  of  Figure  2)  varies  with  the  aperture  of  a  projection  of  the  array  on  a  vertical  plane 
oriented  in  the  target  direction.  For  example,  when  the  target  is  broadside  to  the  16-element  linear 
array,  the  ionospheric  clutter  is  also  broadside  with  respect  to  the  array,  and  the  array  has  difficulty 
distinguishing  the  low-elevation  surface  vessel  echo  from  the  high-elevation  ionospheric  echo.  This 
effect  becomes  more  pronounced  as  the  clutter  spread  angle  a  goes  to  zero,  as  shown  in  the  right  panel 
of  Figure  2.  In  contrast,  the  two-dimensional  arrays  have  better  array  gains  for  targets  near  broadside 
due  to  their  finite  aperture  projection  on  the  vertical  plane  containing  the  target  direction.  The  zenith 
array  null  becomes  perfect  when  a  —  0  and  the  array  gains  for  these  configurations  become  infinite. 
Conversely,  as  a  gets  large,  the  clutter  signal  becomes  uncorrelated  noise  for  all  array  configurations, 
and  the  array  gains  approach  the  uncorrelated  case  of  A  —  N  =  12  dB. 
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Surface  Wave  Propagation 


•  Problem:  diffraction 
of  vertical  polarized 
waves  over  lossy 
sphere  [Fock,  JPM 
1945,  p25 6] 

•  Solve  wave  equation 
in  presence  of  lossy 
boundary 
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Monopole  Array 

•  Transmit  Power 

~16  kW  Pulsed,  10%  Duty  Cycle 

•  PRF 

up  to  250  Hz 
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Range-Doppler:  Spread  F  Conditions 


Time-Range:  Evolution  of  Ionosphere 
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Spatial  Processing:  Beamforming 
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Adaptive  Beamforming 

•  Measured  signal 


y  =  xv  +  n 

•  Estimate 

~  H 

x  =  w  y 


Estimation  error 


e  =  E[(x-x)2] 


•  Goal:  minimize  e  under  constraint 

E(x)  =  x 
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Solution 


*  v^R;1 

w  = 


v"R-'v 


Noise  covariance 


R„=E[nnH] 

Output  noise  power 

R  =  \y/7R  w  = 

x  n 


i 


v"R> 


If  input  noise  is  unity,  then  SNR  gain  of  array  is 

A  =  l/R}  =v"R-‘v 

x  n 


Adaptive  Beamforming  (con’d) 


To  determine  A  for  various  array  configurations,  we  must 

•  Determine  angular  spectrum  of  clutter  signals 

•  Calculate  sensor  covariance  R„  for  an  array  configuration 


•  Determine  array  gain  for  A  for  given  target  location 
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•  Backscatter  occurs  when  surface  normal  to  incident  k 

•  PDF  of  slopes  provides  PDF  of  reflection  angles? 


Slightly  Rough  Reflecting  Surface  (con’d) 


•  By  central  limit  theorem,  PDF  for  altitude  of  reflecting 
layer  at  point  x,  is  Gaussian: 


f(z)  = 


1 


2  <7, 


V2 Tt(7z 

PDF  for  slope  z  ’  can  be  derived: 

1 


f(z')  = 


o  2  2  2 

57T  CJZ  <7k 


(2  TZ)V1CJzGk 


where  crk2  is  second  moment  of  wavenumber  PSD 
Typically,  az’  =  2nazak»  1,  so  not  slightly  rough! 


Turbulent  Plasma 


*,v 


T\0 


6  6 


Ray  tracing  can  be  used  to  determine  PDF  of 
reflection”  angles 


Turbulent  Plasma  (con’d) 


•  Populate  ionospheric  profile  with  irregularities 

•  Choose  an  irregularity  PSD  [Hysell,  JGR  1994,  pi 5078] 


4  71 

An(k±,kn)  =  —^ 


2L.A/7 

0  rms 


ak1  1  +  k2 1^  (1  +  kzL\ ) 


2  t2 


1 


k2  =kl+—k2 

a 


Do  ray  tracing,  find  angular  PSD  of  reflected  signals 
Ray  perturbation  theory  [Witte,  GJI  1996,  pi 65]  shows 


°6  =  i^nns  /  n\\Lpathcrk 


10-1 


z  (m) 


Ray  Tracing  Example 


Ray  Tracing  Near  Reflection  Height 


.x  10 


Angle  Histogram 


Angle  (rad) 


TX 


Measurements 


*,v 


T\0 


Measure  “local  wavenumbers”  Kx(x,y )  and  K{x,y). 

'Ax-*# 

Ay 

A 

ox  dy 

O  is  a  field  with  random  amplitude  a  and  phase  6 : 

<S>(xiy)  =  a(x,y)el9{x,y) 


Since  k2  =  (co/c)2  -  kx2  -  ky2  we  need  only  S(kx,kv): 


S(kx,ky)  =  \\E[®\x,y)®(x  +  X,y  +  Y)]e-i' 


ikxX-iky] 


dXdY 


Measurements  (con’d) 

Assume  slowly  varying  condition,  for  all  m,n>  0 


1  dma 


a  dx 


m 


1  dmK 

A 


K  dx 

A 


m 


« 


K 


m 

x 


1  dna 


a  dyn 


1  dnK 


y 


K,  dy- 


« 


K 


n 


y 


a,  Kx,  Ky  have  joint  PDF  f(a2,  Kx,  Kv): 


S(kx,ky )  *  {Inf  J a\x,y)f{a\Kx,Ky)da2 

•  Wavenumber  PSD  is  power-weighted  histogram  of  local 
wavenumbers 

Analogy:  Frequency  spectrum  is  spectrum  of 
instantaneous  frequencies  if  modulation  rate  «  bandwidth 


Measurements  (con’d) 


*,V 


T\D 


Skymap 


JJU 


S(kx,ky) 


2  K 


G 


kj+k^ 

la1 


spread  «  2cr  /  k 


S(ky ) 


Analysis  of  Array  Configurations 


1 6  Linear  8x8  Cross 


Array  Gain  Versus  Target  Azimuth 


90 


Array  Gain  Versus  Clutter  Spread  Angle 


Clutter  Spread  Angle  (deg) 


Linear  Array:  Broadside  Versus  Endfire 


Broadside  Beamforming 


Endfire  Beamforming 
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Doppler  Freq  (Hz) 


Range  358  km 


Problem  of  Mixed-Path  Propagation 


•  Can  exceed  Bragg  line  power  by  40  dB 


Current  Work 


•  2-D  Array  Experiment 


Conclusions 


•  Goal  was  to  determine  array  gain  A  for  various  configurations 


•  Angular  spectrum  of  clutter  is  Gaussian,  spread  ~  1 5  deg 

•  Array  gains  vary  with  size  of  projection  on  vertical 
plane  containing  target  azimuth 


Cross  array  produces  best  overall  array  gain  for  fixed 
number  of  sensors 


Performance  constrained  by  mixed-path  propagation 
or  array  calibration 
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